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Abstract 

We present a method to prove nonlinear instability of solitary waves in dispersive 
models. Two examples are analyzed: we prove the nonlinear long time instability of the 
KdV solitary wave (with respect to periodic transverse perturbations) under a KP-I 
flow and the transverse nonlinear instability of solitary waves for the cubic nonlinear 
Schrodinger equation. 

1 Introduction 

There are many results (both theoretical and numerical) dealing with detecting unstable 
modes of dispersive equations linearized around soliton like structures. However, in most of 
these cases it is not clear whether one has indeed a nonlinear instability for a flow of the full 
nonlinear problem. The main reason is the lack of understanding of the whole spectrum of 
the linearized problem. The goal of this paper is to present a method showing how only a 
partial information about the spectrum of the linearized operator together with a suitable 
nonlinear analysis may indeed give the proof of the nonlinear instability in the presence of an 
unstable mode. Our first example is the nonlinear long time instability of the KdV solitary 
wave (with respect to periodic transverse perturbations) under a KP-I flow. We also prove 
a nonlinear instability result for the cubic nonlinear Schrodinger equation. We believe that 
the method presented here could be useful in the contexts of other dispersive equations. 

Consider the Kortweg- de Vries (KdV) equation 

u t + uu x + u xxx = 0, (1) 

u : M. 2 — > M, which is an asymptotic model, derived from the free surface Euler equation, for 
the propagation of long one-directional small amplitude surface waves. A famous solution of 
(CQ) is the solitary wave solution, given by 



u(t, x) = Q(x — t), Q(x) = 3 sech 2 
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Observe that u(t, x) corresponds to the displacement of the profile Q from left to the right 
with speed one. One also has the solution 

u c (t, x) = cQ(s/c{x — ct)), c>0 (2) 

which corresponds to a solitary wave with a positive speed c. 

A very natural question concerning the relevance of the solution Q(x — t) is its stability 
with respect to small perturbations. It is evident that the usual Lyapounov stability cannot 
hold because of the translation invariance of the problem. More precisely for c close to one 
cQ(yfcx) is close to Q(x) while for t ^> 1 (t ~ |c — the corresponding solutions of the 
KdV equation u(t,x) and u c (t, x) separate from each other at distance independent of the 
smallness of c — 1. However, the solution u c (t, x) remains close to the spatial translates of Q 
and thus orbital stability of Q under the flow of KdV is not excluded. It is known since the 
seminal paper of Benjamin [2] that Q is orbitally stable in the energy space i/ 1 (IR) (we call 
H 1 (M) the energy space since this is the natural space induced by the Hamiltonian structure 
of ([T])). Here is the precise statement. 

Theorem 1 For every e > there exists 5 > such that if the initial data 

u\ t=0 = u e H\R) 

of the KdV equation (QP satisfies 

\\u - Q\\m(R) < ^ 

then the corresponding solution u (which is well defined thanks to [T^j ) satisfies 

supinf \\u(t,-) -Q(- -a)\\ H i {R) < e. 
tent »gk 

Let us notice that the phase space if 1 (R) may be replaced by L 2 (R) (see [18]). 

In [13], Kadomtsev-Petviashvili studied weak transverse perturbation of the KdV flow 
and derived the following two-dimensional models 

u t + uu x + u xxx ± d~ l u yy = . (3) 

Equation ([3]) with sign + is called the KP-II equation while ([3]) with sign — is the KP-I 
equation. Let us observe that in the derivation of the model, the signs vary in front of the u xxx 
term and correspond to different surface tensions. However from mathematical view point 
the study of the models with signs varying in front of u xxx is equivalent to the study of the 
models with signs varying in front of d~ x u yy by the variable change u(t, x, y) i— > — u{— t, x, y). 
The ant i- derivative d~ l is defined on functions which have, in a suitable sens, a zero x mean 
value. 

Let us observe that Q(x — t) is a solution of both equations ([3]). It is conjectured in 
[13] that Q(x — t) is stable under the KP-II flow and unstable under the KP-I flow. Of 
course this conjecture is very vague since one should precise the stability notion and the 
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spatial domain for x, y. In [TJ, all possible unstable modes of the linearized equation are 
described and in particular it is shown that the linearization about Q of the KP-I flow is 
unstable and the linearization of the KP-II flow is spectrally stable. In this paper, we show 
that the spectral instability result of [1] indeed implies the nonlinear instability in the case 
of the KP-I equation for solutions periodic in the y variable. This result is actually not 
new since the equation being completely integrable (having a Lax pair representation), the 
instability can be shown by exhibiting explicit solutions (see Zakharov [24J). Nevertheless, 
we believe that our method inspired from the work of Grenier [9] in fluid mechanics to 
prove that spectral instability implies nonlinear instability which does not use the complete 
integrability is interesting and can be applied to many other dispersive equations. As an 
illustration, we shall also study below a transverse instability of the two-dimensional cubic 
nonlinear Schrodinger equation which is not completely integrable. 

The global well-posedness of the KP-I equation in the setting IxT was recently ob- 
tained by Ionescu-Kenig [11] in a space which moreover contains the solitary wave Q and 
hence, we state our result in the context of Ionescu-Kenig's theorem. In general it is difficult 
to get nonlinear instability results in natural energy norms like L 2 or H 1 for conservative 
equations due to the presence of strong nonlinearities. Here we shall use the general setting 
developed by Grenier in [SJ in the context of the Euler equation which relies on the possi- 
bility of constructing an high order approximate solution more accurate that the only linear 
approximation. For other methods, we refer to [8], [TO] . 

We consider thus the KP-I equation 

u t + uu x + u xxx - d~ x u yy = (4) 

for (x,y) G M x Tl where Tl is the flat torus M/27rLZ. As mentioned above, a special 
solution of this equation is given by the KdV soliton Q(x — t). Since we are interested in the 
stability of the soliton for (j3J), it is more convenient to go into a moving frame i.e. to change 
x into x — t and to study the equation 

u t - u x + uu x + u xxx - d~ x u yy = 0, (x, y) E R x T L (5) 

so that Q(x) is now a stationary solution of Note that we can always change space and 
time scales to reduce the study of the stability of u c , given by fl2]) to the study of the stability 
of Q for (j3J). Nevertheless, since we are in a bounded domain in y, the scaling changes the 
size of the domain, this is why we keep the parameter L in our study. 

As established in [IT] , the Cauchy problem for (j3J) or equivalently (jSJ) is globally well- 
posed for data in the space Z 2 {R xT L ) defined by 



x T, 



u, \\u(t,k)(i+ iei 2 + ifc/£i 2 )ii L2(RxZ) < +oo}, 



where w(£, k) is the Fourier transform of u : 

™oo p2ttL 



fc ) = 2ttL / / e l ( x ^ +V ^u(x,y)dydx. 



-oo </0 
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If u G Z 2 , this means that u, u x , u xx and d~ l u y , d~ 2 u yy are in L 2 , where d' 1 is defined in 
the natural way via the Fourier transform for functions u G L 2 such that £ _1 w(£, k) G L 2 . 
Moreover, the propagation of H s regularity holds: if u G H s R Z 2 for s > 7, then u(t) G 
n Z 2 for every t > 0. Note that since Q does not depend of y, we have Q & Z 2 . The first 
goal of this paper is to prove the following orbital instability result. 

Theorem 2 Assuming that L > A/\f?>, then for every s > 2, there exists rj > such that 
for every 5 > 0, there exists Uq G Z 2 R H s and a time T 5 ~ | log<5| such that 

— Q\\h s (RxJ l ) < 8 

and the solution u 5 of ([5]) with initial value u 5 satisfies 

inf \\u 5 {T s , ■) - Q(- - a)|| L 2 (RxTi) > rj. 

Remark 3 Ifu(t,x,y) is a solution of the KP equation (Tj|) then so is u\ defined by 

u\(t, x, y) = \ 2 u(X 3 t, \x, \ 2 y). 

Thus in the context of ^ solutions of period one in y transform into solutions of period A~ 2 
and solitary waves of speed c transform into solitary waves of speed \ 2 c. Consequently Theo- 
rem\E implies that if we fix period one perturbations in y then one needs to consider solitary 
waves of sufficiently large speed to get the instability. Let us also remark that the restriction 
L > 4/v^3 in Theorem^ is imposed by the spectral considerations of and is needed for 
the existence of unstable eigenmodes. It would be interesting to decide what happens for 
L < 4/v / 3 (or equivalently for small speed solitary waves for period one perturbations) . 

Remark 4 Let us recall that a three dimensional analogue of ^ 

u t + uu x + u xxx - d~ l (u yy + u zz ) = 0, (x,y,z)eM 3 (6) 

has solutions blowing up in finite time (see fH\j and also [20]) and thus for the three dimen- 
sional versions of the KP-I equation a stronger form of the instability appears. It is however 
an open problem to prove the existence of blow-up solutions for with u periodic in y, z. 

Let us outline the main steps of the proof Theorem [2j First, we need to use the result of [1] 
concerning the existence of unstable eigenmodes for the linearized about Q operator. Next, 
following the idea of Grenier [9], we perform the construction of an approximate solution. 
The approximate solution is defined iteratively. At the first step we put the unstable eigen- 
mode. At each further step, we get linear problems with source terms involving the previous 
iterates (the procedure is closely related to the Picard iteration). We need to control pre- 
cisely the eventual growth in time of each iterate. By applying a Laplace transform, we 
reduce the matters to showing estimates on a resolvent equation which are uniform on some 
straight line A = 7 + ir, r G K. For bounded frequencies (i.e. \t\ bounded), a classical ODE 
argument combined with the absence of unstable modes coming from [lj suffices to get the 
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needed bound. The main difficulty is to get uniform resolvent estimates for large r. They 
will result from conservation (or almost conservation) laws. We finally perform an energy 
estimate to the nonlinear problem to show that the constructed approximate solution is in- 
deed close to the actual solution for suitable time scales. This in turn implies the nonlinear 
instability claimed in Theorem [2j 

The second example that we consider in this paper is the two-dimensional Nonlinear 
Schrodinger equation (NLS) 

iv t + A X)V v + \v\ 2 v = 0, (7) 

where v is a complex valued function. A famous solution of this equation is the solitary wave 
Q(x)e zt with Q given by 

Q(x) = . 

ch{x) 

This solitary wave is orbitally stable when submitted to one-dimensional perturbations i.e. 
perturbations which depend on x only (see [6]). Here orbital stability means that 

Ve>0, 3 5>0 : \\v(0,-)-Q\\ H x m <5 inf \\v(t, ■) - e*V*Q(- - a)\\ mm < e . 



We shall prove that, similarly to the KdV soliton as a solution of the KP-I equation, this 
stationary solution of which is orbitally stable when submitted to one- dimensional per- 
turbation is nonlinearly unstable when it is submitted to two-dimensional perturbation. As 
previously, it is more convenient to set v = e lt u and to study the equation 

iu t + Au — u + \u\ 2 u = 0, (8) 

for (x, y) G R x Tl- A stationary solution of this equation is now given by the ground state 
Q(x). Since the solitary waves modelled on Q(x) for ([7j) are given by 



u\(t, x) = XQ(Xx)e 



i\ 2 t 



we can always reduce by scaling the study of the stability of u\ to the study of the stability 
of Q in (jSJ), but it is again important to keep L as a parameter. Here is our result. 

Theorem 5 There exists Lq such that for L > Lq, for every s > 2, there exists r\ > such 
that for every 5 > 0, there exists u 5 G H s and a time T s ~ | log (5| such that 

\\uq — Q\\h*(}blxt l ) < <5 
and the solution u s of (JBJ) with initial value Uq belongs to C([0, T 5 ], H s ) and satisfies 

inf \\u\T 5 , •) - Q(. - a)e^|| L2(RxTL) > V . 



Remark 6 Let us observe that the cubic two-dimensional NLS is not known to be integrable 
(in the sense of Lax pairs representation) and thus it is hard to expect that the instability 
result presented in Theorem^ can be displayed by an explicit family of solutions in the spirit 

ofM- 
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Remark 7 It is likely that the method presented here may be applied to the following two 
dimensional perturbation of the Boussinesq equation 

u u + (u xx + u 2 - u) xx ±u yy = . (9) 

The stability, for suitable values of the propagation speed, with respect to one- dimensional 
perturbations of the soliton of the Boussinesq equation is obtained in J^j. The analysis for 
an unstable mode in 2D in the context of (Tj|) is essentially the same as the corresponding 
analysis for the KdV soliton as a solution of KP (see fBj]). One thus may perform the analysis 
of |iy (see also the Appendix of this paper) combined with the nonlinear analysis of this paper 
to get statements in the spirit of Theorems^ [3]/br equation |5j). 

The assumption L > L in Theorem is used to get the spectral instability of the solitary 
wave. A difference with Theorem [2] is that for the two-dimensional Schrodinger equation in 
K x T a global existence result of large data strong solutions is not known so that Theorem [5] 
contains the fact that our unstable solution u 5 remains well-defined on a sufficiently long 
time scale. In fact, small data global existence for (181) , posed on K x T, is obtained in [21 J . 
For general large data we may not have the global existence for (jSJ), posed on R x T, since 
one may localize the well-known explicit blow-up solution for the cubic NLS on M 2 (see [5] 
for details on this argument). 

The rest of the paper is organized as follows. In the next section, we give a detailed proof 
of all the steps of the proof of Theorem [2j Then we give a less detailed proof of Theorem [5] 
since the method is the same. Finally, the appendix is devoted to the linear instability results. 



Acknowledgement. We are grateful to Robert Pego for pointing out to us the refer 
ence [23]. A previous version of this text, before we were aware of the Zakharov work [2 
was posted to the arxiv of preprints on December 2006. 



2 Proof of Theorem [2 



2.1 Existence of a most unstable eigenmode 

The linearized equation about the soliton Q reads 

Ut + Au = 0, Au = -u x + (Qu) x + u xxx - d~ l d yy u, (x, y) e R x T L . (10) 

This last linear equation can be solved, for instance by a classical energy method, for initial 
data in H s such that its anti-derivative exists. The main result of [1] is the characterization 
of all the unstable eigenmodes associated to A. An unstable eigenmode is a solution of (fTUl) 
under the form 

i iky 

(p^ k (t } x } y) = e a e~V(x), 
with Re a > 0, V G L 2 (R). The result of [I] adapted to our framework reads : 
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Theorem 8 (pp) There exists unstable eigenmodes if and only if L > A^. Moreover, for 
an unstable eigenmode, a and k G Z are parametrized by 



/or 

2<7 = /i0x-l)(2-/i), fe = X_^( 2 -^), /i6(l,2) (11) 
and there exists g G H°°(R) such that 

V = g xx . (12) 

For the sake of completeness, we recall the main steps of the proof of this result in the 
Appendix. 

Note that for fi G (0,2), /i(2 — /x) G (0,1) hence one can find an integer such that 
k = -^|^ /i(2 — /i) only if > 1. Moreover, for L fixed, there is only a finite number of 
k which verify this property, this allows us to choose <7o and k such that <f ao ,k is the most 
unstable eigenmode i.e. 

a = sup | a, (a,k) verifying ( flTT) j 
and fco is the corresponding integer such that (iTTj) holds with (a, A;) = (<To, fco). Let us define 

u°(t, x, y) = <pao,ko(t, x, y) + Tp^i(t, x, y) . 

To prove Theorem [2], we shall use Q + 5u° as an initial data for (J5j). As remarked before, 
we have Q G Z 2 D -ff 5 for every s, but thanks to (Tl2|) in Theorem [HI we also have that 
u° & Z 2 H H s consequently, thanks to the result of [11] there is a unique global solution u 5 
of ((SJ) in Z 2 fl with initial value Q + <5w°. So the only problem that remains is to estimate 
from below 

inf \\u s (T s ,.)-Q(--a)\\ L * {ML) . 

(i £ M. 

Towards this, we shall use the method of [9] which relies on the construction of an high order 
unstable solution. This is the aim of the next section. 

2.2 Construction of an high order unstable approximate solution 

Let us set v = u s — Q, then v solves 

v t + Av = —vv x . (13) 

We define as the space : 

V£={u,u= u 3 {x)e %3 \ u 3 G # S (R)} 

je-^z, \jL/k \<K 

and we define a norm on V% by 

\u\ V ° = SUp \Uj\ s 

j 



where | • | g is the standard H S (M.) norm. Let us notice that u° is such that u° G V{ for all 
s£N. Following the strategy of [9], for s ^> 1, we look for an high order solution under the 
form : 



M 

ap ' 



5 { u ° + Yl 5kyk ) > uk e v k+ k ( 14 ) 



u 

k=l 



such that Ui t=Q = and M > 1 is to be fixed later. Once the value of M is fixed, then we 
fix the integer s so that s > M. By plugging the expansion in ffl3l) . cancelling the terms 
involving 6 k , 1 < k < M + 1, we choose u k so that u k solves the problem 

8 t u k + Au k = -i( u a=o = - ( 15 ) 

j+l=k-l x 

The main point in the analysis of u ap is the following estimate. 

Proposition 9 Let u k the solution of (jTo"|) . if s — k >1, we have the estimate: 

\u k (t)\ v s- k < C k s e {k+l)uo \ Vt > 0. (16) 

v k+l ' 

The proof of the proposition will follow easily by induction from the following theorem. 
Theorem 10 Consider the solution u of the linear problem 

d t u + Au = F x , u/ t =Q = (17) 
with a source term F 6 V^ +1 with 

\F{t)\ vk+ ,<C F K ^\ 7 >2a (18) 
then u belongs to V£ and satisfies the estimate 

\u{t)\ Vk <C K>s e^, Vt>0. (19) 



We first observe that under our hypothesis on F the solution of ffTUl) is well-defined and the 
only point is to prove the quantitative bound f[T9l . The estimate (|19l) . relies on the fact that 
on Vk, the real part of the spectrum of the operator —A is bounded by o~q. Nevertheless 
for such a dispersive operator, there is no general theory to convert an information on the 
position of the spectrum into an estimate on the semi-group like it is the case for example 
for sectorial operators. To get the result, we need to estimate the resolvent of —A on V^. At 
first, we can perform some reductions on the problem. Indeed, since F has a finite number 
of Fourier modes, we can expand u in Fourier modes and hence we only need to study the 
problem 

d t v + Ajv = d x Fj(t, x), v /t=0 = (20) 

where 

AjV = -v x + (Qv) x + v xxx + fd~ l v } (21) 
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j G ^Z, \j L/k \ < K and v(t, x) = Uj(t,x), and to establish that v satisfies 

=7* 



\v{t)\ s < C jiS e^ 
under the assumption 

<C hS e^. (22) 

In what follows, we fix 7 such that a < 7o < 7 an d we shall use the Laplace transform. 
For T > 0, we first introduce G such that 

G = 0, t < 0, G = 0, t > T, G = Fj, t g [0, T] 

and we notice that the solution of 

fyu + Ajv = G x , v/ t =o = 

coincides with v on [0, T] so that it is sufficient to study 5. Next, we set 

w(t, x) = Cv{^q + ir), H( y T,x) = CG('jo + iT), (r, x) G M 2 

where £ stands for the Laplace transform in time : 

Cf(l Q + ir)= [ e-^- lTt f(t)l t>0 dt. 

We get that w solves the resolvent equation 

(70 + ir)w + AjW = H x . (23) 
In the sequel, for complex valued functions depending on x, we define 

(f,g)= [ f(x)g(x)dx, |/| 2 = ||/||ia(R) = (/,/), I/I^II/IIW)= E i^/l 2 - 

^ M 0<m<s 

Towards the proof of Theorem [10], we first need to study (!23j) . Our main estimate on the 
resolvent will be 

Theorem 11 (Resolvent Estimates) Let s > 1 be an integer. Let w(t) be the solution 
of ( |23l) forj, |j | ^ koK/L, then there exists C(s, 70, K) > such that for every r, we have 
the estimate 

\w(t)\1<C(s, 1o ,K)\H(t)\1 +1 . (24) 
2.2.1 Proof of Theorem HU 

We shall split the proof in various lemmas. To estimate w, we shall deal differently with 
large and bounded frequencies. 

Lemma 12 There exists M > (which depends on K) and C{s ) r jQ 1 K) such that for \r\ > 
M , we have the estimate 

\w(t)\1<C(s, 1o ,K)\H(t)\ 2 s+1 . (25) 
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2.2.2 Proof of Lemma [12] 



We first prove (1251) for s = 1. Note that the equation (1231) can be rewritten as 

( 7o + ir)w - {Cw) x + fd-'w = H x (26) 

where C is defined by 

Cw = —w xx — Qw + w. 

Note that £ is a self-adjoint operator in L 2 which is very useful in the proof of the stability 
of the soliton for the KdV equation. Since it is self-adjoint, the spectrum is real. Moreover, 
since Q goes to zero exponentially fast, the essential spectrum of C is in [1, +oo). For A < 1 
there are only eigenvalues of finite multiplicity. Finally by Sturm-Liouville theory, since Q x 
is in the kernel of C and has only one zero, we get that C has only one negative eigenvalue. 
Moreover, is a simple eigenvalue. Consequently we can define an orthogonal decomposition: 

w = a(r)<p_i + P(r)(p + w ± (27) 

where 

C<P-i = W-i, /i < 0, Cfo = 0, (Cw±,w±) > c |w_l| 2 , c > 0. (28) 

Note that the eigenvectors and (fo are smooth. The important role of £ is due to the 
following conservation law 

7o (fw, Cw) + f\d- l w\ 2 ) = Re Cw) + f(H, d^w)) (29) 

which can be checked by a straightforward computation. Consequently, we can use fl27|) . 
(J2"8"j) and integrate by parts the right-hand side to get 

7o(/ia(r)|^_ 1 | 2 + co\w ± \ 2 + f^wf) < C\H\ 2 \w\ x + f\H\\d x l w\ . 

Therefore, using the inequality 

ab<ea 2 + -^b 2 , Ve > 0, V(a,6)el 2 , (30) 

with e small enough, we can incorporate Ic?" 1 ^! in the left hand-side and arrive at 

\w±\ 2 + f\d~ l w\ 2 < C(\a\ 2 + \H\ 2 + \H\ 2 \w\^j. (31) 

In what follows C is a large number which may change from lines to lines and depend on 7 
and K but not on r. The next step is to estimate a and (3. We use the decomposition (J27j) 
and take the scalar product of ( 1261) with <£>_i and with (p respectively to get 

(70 + ir)a = -(w, Cd x (v-i)) - f(d~ l w, yj-i) + (H x , <p_i) 
(7o + ir)P = -{w, Cd x ((p )) - f(d£ x w, (p ) + (H x , ip ) 
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and hence, we can take the modulus and add the two identities to get 

( 7o + |r|)(M + \/3\) < C(\a\ + \/3\ + \w ± \ +f\d~ 1 w\ + {H^. 

Next, we multiply by \a\ + and use fl30|) to get, 

(70 + \r\ - C)(\a\ 2 + \(3\ 2 ) < C (\w ± \ 2 + f\d~ l w\ 2 + (32) 

Note that this last estimate is a good estimate when r is large. Next, we can consider B (l3TT) 
+ (1321) with B a large number to be chosen to get 

(B - C)(\w ± \ 2 + (Bf - Cf)\d^w\ 2 ) + ( 7o + \r\-C- BC){\a\ 2 + \[3\ 2 ) 

< CB^H^lwlt + lHlfj. 

Consequently, we can first choose B sufficiently large (such that B > C, and B > Cj 2 ) and 
then consider r sufficiently large (for example |r| > 2(C + BC)) to get the estimate 

H 2 + j 2 |d*M 2 <^(l#|2Hi + l#li), \r\>M. (33) 

To conclude we just need to estimate It suffices to look again at Indeed, we can 

use that (w,Cw) = \w x \ 2 - 0(l)\w\ 2 in (J22D to get 

Kl 2 +j 2 \d~ 1 w\ 2 < C(\w\ 2 + \H\ 2 + |fT| 2 |u;|i). (34) 

Consequently, the combination of a sufficiently large constant times (1331) and (1511) gives 

HI +f\d~ 1 w\ 2 < C(\H\ 2 \iv\i + \H\l\ \t\ > M 

and hence by using the inequality fl3"U|) . we get 

IwH+fld-'wf^ClHg (35) 

This proves (|25|) for s — 1. Note that moreover f|35]) gives a control of J 2 1 <9^" 1 w | 2 which is 
interesting when j ^ 0. 

To estimate higher order derivatives, we shall use higher order approximate conservation 
laws for the linearized KdV equation. Namely, we define a self-adjoint operator 

C s+1 w = d 2 x s+2 w + r s+1 {x)d 2 x s w 

where r s+ i is real valued and will be chosen in order that the following cancellation property 
occurs : 

Re (iCw) x , C s+1 w^) =0{l)\w\ 2 . (36) 
11 



By making repeated integration by parts, we easily establish that 



Re (d 2 x s+2 w, w xxx ) 
Re (d 2 x s+2 w, Qw x ) 

Re (d 2 x s+2 w, Q x w) 
Re (r s+1 dl s w, d xxx w) 



l) s Re (d s x +2 w, d x d s x +2 w) = 0, 

(d£+V Qd s x +2 w) + (s + l)(d s x +1 w, Q x 3 s x +1 w) + 0(1) \w\ 



-l) s+1 Re 
-l) s+1 Re 
-l) s+1 Re 
-l) s - x Re 
-l) a_1 Re 



[s + -)(d x +1 w,Q x d s x +1 w) + 0(l)\ 
(d s x +1 w,Q x d x +1 w)+0(l)\w\ 2 ), 



w\ 



d x +1 w,d s x 1 (r s+1 w xxx 



(d s x +1 w,r s+1 d s x +2 w) + (s- l)(d s x +1 w,d x r s+1 d x +1 



w , 



+ o(i)K 

3 



ly-'Re ((, - ^m +1 w, d x r s+1 d x +1 w) + 0(l)\w\l) 



and that all the other terms which appear in the product Re ((Cw) x , C s+ iw) are 0(1) Iw^. 
Consequently, we get 



Re(-(£w) x ,£ s+1 w) = (-iy^(( s + -)(d s x +1 w,Q x d x 

+ (s-lm +1 w,d x r s+1 d x +1 w)) 



+ 0(l)\w 



0(l)\w 



with the choice 



r s +i 



Q. 



Note that s is an integer so that r s+1 is always well-defined. 

Finally, we can take the scalar product of (1261) by (— l) s+1 £ s+ iw and then take the real 
part to get thanks to the above cancellation property 

lo\d s x +1 w\ 2 < C(\w\ 2 s + f\d x l w\ 2 + \H\ s+2 \d s x +1 w\ + |#| s+1 Ms) 

since Re {d~ x w, dl s+2 w) = 0. We finally obtain 

lW<c(M!+jW + l^t 3 ) 

thanks to the inequality fl30l) and hence we get (1251) by induction and the control of j 2 \d~ 1 w\ 2 
given by fl35l) . 



Next, we need to estimate w for |r| < M. This is the aim of the following lemma. 
Lemma 13 For Irl < M , we have the estimate 



\w(t)\ 2 <C(s, 1o ,K,M)\H(t)\': 



+1- 



(37) 
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2.2.3 Proof of Lemma [H 

Note that here we actually give a proof of the fact that if A is not an eigenvalue then A is 
not in the spectrum. To prove fl37|) . we need to treat differently the cases j = and j ^ 0. 
Let us start with the case j ^ 0. In this case, we take the derivative of fl26|) to get 

(70 + ir)w x - (Cw) xx + fw = H xx (38) 

and we introduce V = (w, w x , w xx , w xxx Y G C 4 and H = (0, 0, 0, H xx ) 1 to rewrite the problem 

as 

V x = A(q,x)V + B. (39) 

where A is a 4 x 4 matrix that one may easily find from the equation (|38il and the parameter 
<? = (7o + ii~,j 2 ) is in the compact set /C defined by 

IC = {(70 + ir, b), \t\ < M, k 2 /L 2 < \b\ < K 2 k 2 /L 2 }. 

Let us denote by T(q, x, x') the fundamental solution of 14 = AV i.e. the solution such that 
T(q,x',x') = J4. Next, since Q(x) tends to zero exponentially fast when x — > ±00, there 
exists a matrix A 0O (g) such that 

A(q,x) - Aoo(g) = 0(e~ lxl ), x -> ±00. 

Moreover the eigenvalues of Aoo are the roots of the polynomial P defined in fl63l) below 
and hence are not purely imaginary. By classical arguments of ODE (namely the roughness 
of exponential dichotomy, see [7] for example), the equation V x = AV has an exponential 
dichotomy on M + and M_, i.e., there exists projections P + (q,x), P~(q,x) which are smooth 
in the parameter with the invariance property 

T(g, x, x')P ± {q 1 x') = P ± (?, x)T(q, x, x') (40) 

and such that there exists C and a > such that for every U G C 4 , and q G /C, we have 

\T{q,x,x')P + (q,x')U\ < Ce"*^-^ \P + (q,x')U\, x > x' > 0, 
\T(q,x,x')(I - P + (q,x'))U\ < Ce^ x ~ x,) \{I - P+(q,x'))U\, < x < x', 
\T(q,x,x')p-(q,x')U\ < Ce a{x - X ' ] \p-(q,x')U\, x < x < 0, 
\T{q,x,x')(I - p-{q,x'))U\ <Ce~ a{ - x -^ \{I - P~(q,x'))U\, > x > x' . 

In particular, note that a solution T(q, x,0)V° is decaying when x tend to ±00 if and only 
if V° belongs to lZ(P ± (q 1 0)). Since by the analysis of pQ recalled in section |4~T1 there is no 
eigenvalue of Aj (see (1211) for the definition of Aj) for q G /C, we have no non trivial solution 
decaying in both sides and hence we have 

TZ(P + (q,0))nlZ(p-(q,0)) = {0}. (41) 
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Let us choose bases (rf,rf) of TZ(P ± (q, 0)) which depends on the parameters in a smooth 
way (see [13] for example) then we can define 

M(q) = (r+,r+, 77,7^) 

and we note that M(q) is invertible for q G /C because of ( 1411) . This allows us to define a 
new projection P(q) by 

P(q) = M(q)( I * o)M(g)- 1 

and next 

P(q,x)=T(q,x, O)P(q). 

The main interest of these definitions is that we have TZ(P(q)) = TZ(P + (q,0)) and 1Z(I — 
P(q)) = 7Z(P~(q, 0)). Therefore thanks to (14"U|) . we have for every x that lZ(P(q,x)) = 
TZ(P + (q,x)) and similarly that 

K(I-P(q,x))=Tl(P-(q,x)). 

Consequently we have the estimates 

\T(q, x, x')P(q, x') \ < Ce"^ 1 "^ ,x,x'eR,x> x', Vg G /C, (42) 
\T{q,x,x'){I - P{q,x'))\ < Ce^"^ , x, x' e R, a; < a/, Vg G /C. (43) 

By using this property, the unique bounded solution of (1391) reads by Duhamel formula 
T(q, x, x')P{q, x')U(x') dx - / T(q, x, x')(I - P(q, x'))H(x') dx 

■OO J X 

and hence, we get thanks to f|42l . (|43|) that 

which yields by standard convolution estimates 

\V\ < C\H\. 

The estimates of high order derivatives is very easy, it suffices to write 

d s x +1 V = Ad s x V+ [d s x , A]V + d x U, 
and to write Duhamel formula considering [d^, A]V as part of the source term. 

It remains the case j = 0. In this case, we do not take the derivative of (|26|) . we directly 
define W = (w,w x , w xx ) and we rewrite (126|) under the form 

W x = M(X,x)W + M. 

Then the proof of the estimate follows the same line, we find that has no eigenvalue on 
the imaginary axis. This yields that there is an exponential dichotomy on M + and R_ for 
this system. Next since, the spectrum of the linearized KdV equation about the soliton is on 
the imaginary axis, we get that the system has an exponential dichotomy on the real line. 
We do not detail more since the proof is similar to the previous case. 
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2.2.4 End of the proof of Theorem [TT] 

To get (12"4"1) , it suffices to combine Lemma [12] and Lemma [J3j 



2.2.5 End of the proof of Theorem EH 

By using Theorem [11] and Bessel-Parseval identity, we get that for every T > 0, 



T p+oo 

e- 2 ^ 0t \v(t)\ 2 s dt< / e - 270 '|S(t)|^t = / \w(r)\ 2 s dr 

</0 



C / \H(r)\ 2 s+1 dr= [ T e- 2 ^\F,(t)\ 2 s+1 dt 
Jr Jo 

and finally thanks to f[T8|) . we get 

e~ 2lot \v(t)\ 2 s dt<C e 2 ^ ,o)t dt<Ce 2{ ^ ,o)T (44) 



< 



since 70 was fixed such that 7 > 70- To finish the proof, we notice that the energy estimate 
for the equation f[20|) gives 

Consequently, we can multiply the last estimate by e~ 27oi and use f[22|) to get 

^(e- 2 ^Kt)| 2 ) <c(e- 2 ^\v(t)\ 2 s + e 2 ^ o)t ). 



Next, we integrate in time and use (1441) and again the fact that 7 > 70, this yields 

e~ 2l0t \v(t)\ 2 < Ce 2 ^ o)t . 
This ends the proof of Theorem [10] . 

2.2.6 Proof of Proposition O 

By induction, it suffices to use Theorem [TO] and the fact that H S (R) is an algebra for s > 1. 

2.3 Nonlinear instability: end of the proof of Theorem [2] 

Of course, we only need to prove the statement for 5 small enough. Let us define w by 
setting v = u ap + w, where u ap is defined by f[T4|) . Therefore we have that the solution u s 
may be decomposed as follows 

u 5 = Q + u ap + w . 



If we set 



F = {d t + A)u ap + u ap u 
15 



ap 

x 1 



8 _ Mk/5) 



where A is defined in (fTUj) , then thanks to Proposition [9J 

||F(tr)l|L 2( MxT L )<CM^ M+2 e( M+2 ^. 
We have that u> solves the problem 

(d t + A)w + d x (u ap w) + ww x + F = , w /t=0 = 0. (45) 
We now estimate the solution of (1151) . Using that 

7mxT l 

multiplying (|45[) by w and integrating 1 x T^, we get after several integrations by parts 

j t \Mt, -)\\b < (Wh~ + \\d x u ap (t, + l) \\v,{t, + \\F(t, ■ (46) 
Observe that 

M 

RMV)IU-(ttxT £ ) < E Cfc .^ +le( " +1)ffo1 ' 

fc=0 

Next, we set 

rpO — 

where k g]0, 1] is small enough to be chosen after the several restrictions we will impose in 
the next lines. The number T s represents the time when the instability occurs. Coming back 
to (jl6j) . we observe that there exists a constant A MjS depending on s and M but independent 
of k and t and an absolute constant C (C is essentially ||Q'||l°°) such that for < t < T 5 , 

j t \Mt, <(C + KA M , s )\\w(t, -)|| 2 L2 + C M ,^<"+V<" + ^* . 

Therefore 

^ e -(-A«, s +^|| w(t) .)||2\ < c , Mis5 2(M +2 ) e 2(M+2) 0i) t-«A M „ t -Oi ^ t G [Q)T 5] (47) 

Now we choose M large enough and k small enough so that 

2(M + 2)ct - kA a/jS - C > 0. 

At this place we fix the value of M (and of s, for exemple s — M + 1) while we will make 
two more restrictions on k. Since w vanishes for t = an integration of (1471) yields 

IN*, OII^CRxT*) < CM, s 5 M+2 e (M+2)<TO ' , t G [0,T 5 ] . 

Therefore 

lh(T 5 ,-)ll^(Kx^) <C M)S k m+2 . (48) 
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Let us denote by II the projection on the nonzero modes in y i.e. 



2nL 



(Uv){x,y) = v(x,y)-^-J v(x,y)dy. 

Then for every a G R one has U(Q(x — a)) = 0. On the other hand the first term of u ap 
satisfies I1(m ) = u° and therefore 

M M 
k=l k=l 

where c s is the H S (R x T/,) norm of u°. Therefore for k small enough one has 

\\U(u^(T s ,.))\\ L2{ML) > C -f. (49) 
Using PHD and ( 1491) . we may write that for every a G R, 
\\u s (T s , ■) - Q(- - a)|| L2 > \\U(u 5 (T s , ■) - Q(- - a))\\ L , 

= linger 5 , •) - Q(.))\\ L2 = \\u(u^(t s , •) + w(t 5 , 0)11* 



> ^- pht 5 , .))||* >f -HrV)b > ^-c My , 



Af+2 



A final restriction on k may insure that the right hand-side of the last inequality is bounded 
from below by a fixed positive constant r] depending only on s (in particular 7] is independent 
of 5). This completes the proof of Theorem [2j 

Remark 14 Let us observe that the analysis in the proof of Theorem\^is quite different from 
the high frequency instabilities studied in [16]. In fTb^ . the approximated solution is a high 
frequency linear wave with modified speed, perturbed by a low frequency wave. In Theorem® 
the approximated solution is a low frequency object modelled on the profile u°. 



3 Proof of Theorem [5 

The proof follows exactly the same lines as the proof on Theorem |2] and thus we shall only 
sketch it. We again look for u s under the form u 5 = Q + u ap + w. At first, we need to find a 
most unstable eigenmode for the linearized equation to begin the construction of u ap . The 
linearized equation about Q reads 

iut + Au = 0, Au = Au - u + 2uQ 2 + uQ 2 . 

It is more convenient to introduce U = (Rett, Imit)* and to rewrite the equation as the 
system : 

U t +( £ ° + "o " ) U = 0, (50) 
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C u = —Au + u — Q 2 u, C + u = —Au + u — 3Q 2 u. 
We seek unstable eigenmodes under the form 

i iky — i — iky 

$ a>k (t,x,y) = e at e— V(x) + e at e— y(x), Recx>0, (51) 
where V(x) G C 2 so that we have to solve 

+ g ~ L \~ % ) v -° < 52 > 

where 

L~u = —u xx + u — Q 2 u, L + u = —u xx + u — 3Q 2 u. 

We set £ = j and we look for nontrivial solutions of fl52|) with Re<r > for e > 0. The first 
result we shall use is that 

Lemma 15 For e > 0, there is at most one unstable eigenmode and there exists Eq such 
that for < e < Eo, there is exactly one unstable eigenmode. 

In the reference [12], it is claimed that the result of this lemma is due to Zakharov-Rubenchik. 
Unfortunately, we were not able to find a copy of the paper by Zakharov-Rubenchik as this 
paper is quoted in [12] . We give a proof of this Lemma in the appendix. 

Now, thanks to Lemma [TSl for k — 1 and L sufficiently large there exists an unstable 
eigenmode. We now consider L as fixed. For every k, we have by Lemma [T5l that there exists 
at most one o~(k) such that Rea(k) > and (|52l) has a solution in L 2 (IR; C 2 ) with a = a(k). 
Moreover we can easily get that the solutions of (|52|) satisfy the conservation law 

Rea((L+^, V x ) + (L~V 2 , V 2 ) + ^\V\ 2 ) = . 

Therefore for large k (depending only on Q) there is no nontrivial solution of fl52l) with 
Re a > 0. Consequently, we can choose an eigenmode <& a ^ under the form ( ISTI) such that 

Rea = supjRe a(k) j := a 

and we set u° = ($cr,fe)i + i(&a,k)2- Observe that thanks to (|52l) we have {id t + A)u° = 0. 
The next step towards the proof of Theorem [5] is the construction on an high order unstable 
solution. We use the same method as previously, we use the same spaces V£ and we build 
an approximate solution under the form (TH1) . For 1 < k < M + 1, we need to solve 

id t u k + Au k = - {ZQu 3 u l + Qu j u 1 )- u j u l u m , (u k ) /t=0 = (53) 

j+l=k-l j+l+m=k~2 

where the last sum is zero for k = 1. We have the estimates : 
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Proposition 16 Let u k the solution of (|53|) . we have the estimate 

W k {t)\ V s +i <C k J k+1 ^ \ Vt >0. 

Note that here we do not loose regularity at each step because the nonlinear term does not 
involve derivatives. To prove Proposition [16j we need to prove the equivalent of Theorem [TUl 
By using Laplace transform, we can still reduce the problem to the proof of a resolvent 
estimate as in Theorem [TH The proof of the low frequencies estimates rely on the same ODE 
argument and we shall not detail it. We shall just explain how to get the high frequencies 
estimates. As in Lemma [TBI it is more convenient to work on the system form of the problem, 
and thus we consider the equation 

(l + iT)W+( L+ ° + p ~ L ~-^^W = H (54) 

and we want to prove that W(r) satisfies the estimate 

\W(r)\l<C(s, l0 ,K)\H(r)\l (55) 

for 70 > (To, \t\ > M ^> 1 and s > 1. We first give the proof for s = 1. The conservation 
law reads for W = (w±, w 2 ) 

1 2 

lo ({L+w l ,w 1 ) + {L~w 2 ,w 2 ) + —I W\ 2 ) = Re ((# 1; L + Wl ) + (# 2 ,L^ 2 )). (56) 

At this stage, we shall use the description of the spectrum of L recalled in the appendix of 
this paper. We can write 

W2 = aQ + w 2 , (L~W2, w 2 ) > cq\w 2 | 2 . 

Similarly, we can write 

Wi = /3(p-! + 'jQx + w^, (L + w^, w^) > c \wi\ 2 
(ip_ x = Q 2 ). Setting W 1 = (w^w^Y and W F = (a, £,7)' G C 3 , we get from (j55]l 

7o|W^(t)| 2 < cflH^ \W\! + \W F \ 2 y (57) 

Next, we can take the projection of the equation on the finite dimensional subspace generated 
by (0,Q), (Q„0), (<^_!,0) to get 

(70 + |t| - C)\W F \ 2 < C{K)(\W L \ 2 + |#| 2 ). (58) 

As for the KP-I equation, a suitable combination of (156]) . (joTj) . fl58l) with the use of fl30l) 
gives fl55l) for s = 1 for |r| large enough. To get higher order derivatives, we use approximate 
higher order conservation laws. Namely, we choose Lf +l and L~_ x under the form 

L- +l w = d 2 x ^w + r 7+1 d 2 x s w Lf +1 w = d 2 x ^w + r+ +1 «9> 
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such that the cancellation 

-Re (l"w 2 + — W2 ,Lt +1 w^ + Re + _ Wl) L; +1 ™ 2 ) = + |<9* +1 jy|) 

occurs. To perform this cancellation, it suffices to choose rf +1 real valued and such that 

rt + i-r- +1 = 2Q 2 . 
Using this approximate conservation law, we get 

7o|^ +1 W| 2 < c(\W\ 2 e + \W\ S \W\ S+1 + \H\ S+1 \d s x +1 W\ + \H\ S \W\ S ) 



and we conclude thanks to ( 1301) via an induction argument. 

To end the proof of Theorem [51 we seek for a solution of ([8]) under the form u s = 
Q + u ap + w, with w/ t =o = so that w solves the equation 

iw t + Aw + 2\u ap \ 2 w + (u ap fw + AT{u ap , w) + \w\ 2 w = F (59) 

with 

\\F\\h^t l ) < C M ,5 M+2 e (M+2) °°\ 
and the bilinear term satisfies 

\\Af(u ap ,w)\\ H s^ xTL) < C\u ap \ W s^\\w\\ 2 s . 

Since here we do not have a global existence result available, we shall first prove that this 
last equation has a smooth solution w which remains defined on a time scale sufficiently long 
to see the instability. 

A classical existence result for this equation based on Duhamel formula and Sobolev 
embedding gives that there exists a local solution w G C([0,T], H s ) for s > 1. Moreover, we 
can define a maximum time T* such that 

T* = su P {T,Vt e [0,T],\\w(t)\\ H s < I}. 

The H s energy estimate for (159]) gives for t 6 [0, T*) that 

j t \\w{t)f HS < C(l + \u ap \ w ^)\\w\\ H s + CM^^e 2 ^ 2 ^ 

where C is an absolute constant (which depends on Q). Consequently for 

t< Min(T 5 :=^^),n, 

we get 

j t \\w(t)\\% <(C + KA^HIfl- + C M ^ M+2 h 2 ^ M+2 ^ 
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and hence by the choice 

2(M + 2)tr - kA m<8 - C > 0, 

we get that 

\\w(t)\\ Hs(Mt) < C m>s k M+2 , t < Min (T 5 ,T*). (60) 
In particular for k sufficiently small, we get that 

\Ht)\\ H s [ML) < i t< Min (T S ,T*). 

By definition of T*, this proves that T* > T s so that the time of existence of a smooth 
solution is in any case large enough to see an instability. The end of the proof follows the 
same lines as previously, using again the projection IT on nonzero modes in y, we write for 
every a 6 R, 7 6 I, 

|| U *( T *, ■) - e *Q(. - a)\\ L2 > \\U{u s (T 5 , •) - e *Q(. - a))|| L2 = ||n( M ^(T 5 , •) +w(T s , 

> ^ - ||n( w (T 5 , .))||*. > ^ - || W (T 5 , 01k- > ^ - ^ M+2 • 

where we have used ( 160]) in the last inequality. A final restriction of k gives the instability 
result. 

4 Appendix 

4.1 Proof of Theorem [8] 

In order to have the same equations as in [1] , we look for solutions of ( fTUl) under the form 

u{t,x,y) = e^e^U^j 

with [/ e L 2 , Re A > and k 7^ 0. Note that this last condition is natural since for k — 0, 
we cannot find instability since the KdV soliton is stable in the KdV equation. We get for 
U the equation 

AXU Z + 4($ 17) M + £/« M - 4tf M + 3t7 2 C7 = (61) 

where we have set 

9 16P , . 

= -jt (62) 

and $ = 3 sech 2 ^. Since $ and its derivatives tend to zero exponentially fast when z — > ±00, 



the solutions of ( ISTT) have the same behaviour as the solutions of 

AXU Z + U zzzz - AU ZZ + 3 V 2 U = 

when z — > ±00. The characteristic values /x of this linear equation are the roots of the 
polynomial P defined by 

P(jj) = /i 4 - V + 4A/i + 3r/ 2 . (63) 
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Consequently for 77 7^ and 7 = Re A > 0, fx £ iR. Indeed, if fx = i£ G iR, then £ should 
solve 

£ 4 + 4£ 2 + 4A£i + 3r/ 2 = 

which cannot have a real root £ for 77 7^ and ReA 7^ 0. A consequence of this is that the 
number of roots fx of positive real part of P is independent of the parameters. Since the 
limit 77 — > +00 gives 

// = 3*0^77 + 0(1), w 4 = -l 

we finally get that P has two roots of positive real parts and two roots of negative real parts. 
This proves that the solutions of (16T|) either tends to zero or blows-up exponentially fast 
when z — > ±00. Moreover, the stable manifold and the unstable manifold have the same 
dimension 2. Finally, there will be a nontrivial bounded solution of (I6TI) if and only if U 
belongs simultaneously to the stable and the unstable manifold. 

In our case, this condition can be computed explicitly. Indeed, we notice that for 7 > 0, 
77 7^ there is a bounded solution of (16T1) if and only if U — g zz with g bounded which solves 

9zzzz + 4 $ g zz + 4 \g z - Ag zz + 3rf 2 g = (64) 

Note that the asymptotic behaviour of the solutions of this equation is also determined by 
the characteristic values given by the roots of P so that this equation also has stable and 
unstable manifolds of dimension 2. Moreover, if ft is a root of P, then 

9fl ( z ) = (V + 2fi + A - 3fi 2 tanh zj (65) 

is a solution of (161)1 . In particular, if Re fx > 0, then is in the unstable manifold. Moreover, 
when P has two simple roots fii, fi<i of positive real parts, then one can prove (see [1] for 
details) that g^ 2 are linearly independent so that they constitute a basis of the unstable 
manifold. Consequently, any bounded solution of (1611) must be a linear combination of g^, 
and . 

Now, let us define 

C+(n)= lim e-^ z g^ = fx 3 + 2fx + \-3fx 2 . 

z— >+oo 

Then, if C + {fXi) 7^ 0, i = 1, 2, we cannot have nontrivial solutions which tend to zero when 
z — > +00. Consequently, this proves that when the positive real part roots of P are simple, 
then a necessary condition to have bounded solutions of (161)) is that C + (fi) = for some 
root fx of P of positive real part. In the case where fx is a double root, then one can check 
that the same condition holds. Indeed it suffices to take and d^g as a basis of the unstable 
manifold (again, we refer to [T] for details). 

It remains to study the equation C + (fx) = with fx a root of P of positive real part. This 
yields the system of algebraic equation 

P(fx) = 0, fx 3 + 2fx + A - 3/i 2 = 0, (66) 
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with the constraint Re fi > 0. The elimination of A between the two algebraic equations 
gives 

A = -^-l)(//-2), r/ 2 = /i 2 (/i - 2) 2 . (67) 

The analysis of this system gives that there is a solution with Re A > 0, Re/i > 0, if and 
only if given fi G (0, 2), 77 and A are given by 

T) = n(2-fi), A = -fi(ji- l)(//-2). 

Finally, we notice that when C + (/i) = 0, we have 

g^z) = 3/jV z (1 - tanh z) = 0(e- (2 -^) 

and hence lim z _ ++00 ^ = since 2 — \i > 0. This proves that C + (/i) = with /x a root of 
P of positive real part is also a sufficient condition to have a bounded solution on R. This 
ends the proof. 

4.2 Proof of Lemma [15 

Set V(a;) = v(x)y with u, t> real valued functions. Then fl52|) implies that 

L + m + £ 2 u = — av , + e 2 v = au . (68) 

Observe that if (u, v) is a solution of (|6"8~|) corresponding to a complex number a then (w, — u) 
is a solution of (1681) corresponding to —a. The operators L + and L~ have classical self adjoint 
realizations on L 2 (M) and their spectrum are well-known (see e.g. [221 [23]). The operator L + 
has exactly two simple eigenvalues —3 and with corresponding eigenfunctions Q 2 and Q'. 
The continuous spectrum of L + is [1, oof. The operator L~ has only the simple eigenvalue 
with corresponding eigenfunction Q and the continuous spectrum of L~ is [l,oo[. Observe 
that (I68p may be written as 

^M:;)(r I , + ,)(;)-(;)' m 

Thanks to the above discussion on the spectrum of L + and L~, we obtain that 

( L+ + e 2 \ 
V L-+e 2 J 

has at most one negative eigenvalue which should be simple. Therefore, thanks to [T0[ 
Theorem 3.1], there cannot be more than one unstable mode. 

For e < 1, the bifurcation of the eigenvalue zero in the case e = can be explicitly 
computed. Note that zero is an isolated eigenvalue so that we can use perturbation methods 
as in finite dimension (see [13] Theorem 1.8, Chapter 7). In the case e = 0, we have that 
zero is an eigenvalue of multiplicity 4 for the linear map introduced in the left hand-side of 
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(169|) (see [23]). The generalized eigenspace splits into two two dimensional invariant sub- 
spaces corresponding to the eigenvectors (u,v) = (Q',0) and (u,v) = (0,Q) respectively. As 
generalized eigenvectors, we can take \(Q + xQ x , 0) and (0, \xQ) which verify 

\(Q + xQ x )\__( 0\ £ f0 \_fQ a 



/ V Q J V l x Q J V 



Thanks to the analytic dependence in e (see [IS]), we look for a a in ( |68|) of the form a = 
iij\E + oj 2 e 2 H — • with Re (u>i) > which corresponds to an unstable mode. We will see below 
that the invariant subspace corresponding to (u,v) = (Q',0) splits to two one-dimensional 
invariant spaces corresponding to eigenvalues with uoi purely imaginary and, what is of 
importance for our purposes, the invariant subspace corresponding to (u, v) = (0, Q) splits to 
two one-dimensional invariant spaces corresponding to eigenvalues with positive and negative 
u)\. The eigenvector corresponding to a positive ui\ provides the unstable eigenmode. Assume 
that u and v are expanded as 

U = Uq + U\E + U 2 e 2 H , V = V + ViE + v 2 e 2 H . 

Then (uq, Vq) satisfy L + (uq) = L~(vq) = 0. Thus there exist two numbers a® and f3$ such that 
Uq = a Q' and v = (3qQ. Then (ui,Vi) are solutions of L + (ui) = —loi(3qQ, L~(v\) = uia Q'. 
Therefore there exists two numbers ot\ and (3\ such that 

ui{x) = ^(xQ'(x) + Q(x)) + ai Q'(x), v x {x) = -^-( x Q(x)) + (3 x Q(x) . 

Next, {u2-,v<i) are solutions of 

L + (u 2 ) = -a Q' - ut ( - -y^-(xQ) + (3iQ) - u 2 f3 Q Q, 

L- (v 2 ) = -(3 Q + ojx (^-(xQ' + Q) + a x Q') + uj 2 a Q' . (70) 



The first equation of ( 1701) can be solved if the right hand side is is orthogonal to Q' (the 
kernel of L + ). This imposes that either «o = or 

- a Q'(x) -a>i^ — -y^(xQ(x)) + PiQ(x)J - cu 2 p Q(x)^Q'(x)dx = 0, 

which implies that uj\ = —49 2 , where 9 = \\Q'\\l 2 (m.)/\\Q\\l 2 (m.), i-e. u)\ = ±i9. Hence if ao ^ 
we have an eigenmode with purely imaginary uj\. 

The second equation of fITOl) can be solved only if the right hand side is orthogonal to the 
kernel of L~ , i.e. to Q. This imposes that either /5 = or 

PoQ(x) + tu 1 (^^-(xQ'(x) + Q(x)) + otiQ'(x)) + uj 2 a Q'(x)^Q(x)dx = 
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which implies that io\ = 4, i.e. ui\ = ±2. From the above discussion, we have that either 
a = or p = 0. If a 7^ (and thus /3 = 0) we obtain purely imaginary oj\ and have the 
bifurcation of (Q' , 0). These modes are not of interest for us. If (3q ^ (and thus a = 0) we 
indeed have en eigenvalue with positive U\. This mode corresponds to the eigenvector which 
is the bifurcation of (u, v) = (0, Q) to the unstable mode of the form floTl) for the linearized 
about Q cubic NLS equation. 
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